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SUMMARY 


A profile of the ice cover in the southern Beaufort Sea was 
obtained by the nuclear submarine USS GURNARD from 7 to 10 April 1976, 
using a narrow-beam upward-looking sonar. The 1400 km profile 
consisted of three legs, of which the long south-north and east-west 
legs intersected near the Caribou camp of the AIDJEX experiment. 

After initial processing at the Arctic Submarine Laboratory, San Diego, 
the data were passed via the AIDJEX Office to Scott Polar Research 
Institute for statistical analysis. 


The analysis was carried out over contiguous 50 km sections of 
the profile, and yielded the following significant results:- 


1. Probability density function of ice draft. The profile was found to 


consist of three distinct types of ice cover. Type 1, nearest the coast 
and extending from the 100 m isobath northwards to 71 30'N at 144°13'W, 
consisted of heavily ridged ice with a mean draft of 4.2-5.1 m. Type 2 
comprised the major part of the profile (1200 km), and was shown 
statistically to be a homogeneous ice cover with the following distribution 
of ice drafts:-— 0-0.5 m, 1%; 0.5-2 m, 10%; 2-5 m, 72%; over 5 m, 17%; 
mean draft 3.67 m. Type 3 was another short stretch of track at the far 
western end of the profile (151°W to 154°14'W at 72°42'N) which again 
consisted of more heavily ridged ice of mean draft 4.5—-4.6 m. The 
proportion of thin ice was the most variable quantity in the distributions, 
having a range 0.1-3.5% (0-0.5 m draft) and 0.4-12.3% (0-1 m draft). The 
cumulative probability distributions for two 200 km sections crossing the 
AIDJEX camps were calculated for application to the AIDJEX model. 


ae 


undeformed floes contaminated with a few keel bottoms - was found to comprise 
56% of the ice in the type 2 cover. The most probable draft of level ice 
was 2.7-2.8 m, with secondary peaks at 2.0-2.1 m and 3.1-3.2 m. Thin ice 
had a peak at 0.8-0.9 m. 


2. Level ice. Ice of less than 1 in 40 local gradient - i.e. mainly 


3. Keel spacings. The distribution of keel spacings was found to obey a 
negative exponential distribution, with departures at small spacings (a 
deficit due to a "keel shadowing" effect, where shallow keels are concealed 
by neighbouring deep keels) and at large spacings exceeding 300 m (a surfeit 
due to the presence of leads and polynyas). 


4, Keel drafts. The distribution of drafts did not obey the theory of 


Hibler et al (1972) which predicts P(h) exp(-bh?), but rather a simpler 
distribution of form 


P(h) dh = B exp(-bh) dh 


with B, b as parameters. Hitherto only sails have been found to obey 

this relationship, and it is postulated that the cause is the very narrow 
beamwidth of the transducer, which enables it to probe into the fine 

structure of a keel and thus perceive multiple secondary "keels" where a 
wide-beam sounder sees only a single entity. Again types 1 and 3 have heavier 
ridging than type 2. For keels deeper than 9 m, type 2 ice has a mean of 

1.39 keels per km with a mean draft of 11.79 m; types 1 and 3 have up to 

5.1 keels per km with a mean draft of 12.1-12.7 m. Maximum keel draft was 
31.12 m, with only one keel over 30 m in draft. 


5. Leads and polynyas. A lead was defined as a continuous sequence of 
depth points where no point exceeds 1 m in draft. The average spacing of 
leads was 212 m over the whole profile, but almost all leads were less than 
50 m in width. The mean spacing of leads over 50 m wide was 10.3 km, and 
of leads over 500 m wide was 237 kn. 


1. INTRODUCTION 


During the period 7 to 10 April 1976 the nuclear submarine 
USS GURNARD (SSN-662) obtained a sonar profile of length 1400 km 
under the Beaufort Sea ice cover in the vicinity of yl AIDJEX 
(Arctic Ice Dynamics Joint Experiment) main camp. re 1 shows 
the route followed by GURNARD, comprising cieTiuerte fi ws and 
east-west (RPS) legs intersecting near the Caribou camp (C), 
together with a connecting leg QR. GURNARD was equipped rae a 
high-frequency, narrow-beam, upward-looking sonar installed by the 
Arctic Submarine Laboratory, Naval Undersea Center, San Diego. The 
sonar fed its output into a signal processing system that digitised 
the range to the ice underside, subtracted this from the transducer 
depth and thus generated a digital magnetic tape of ice drafts with 
a nominal resolution of 0.03 m. Mechanical limitations within the 
system reduced the absolute accuracy of ice draft to + 0.3 m, with 
a demonstrated standard deviation for smooth ice of 0.09 m. Other 
details of the sonar system and of the submarine's depth and speed 
remain classified. The Arctic Submarine Laboratory did, however, 
supply information on the "surface beam diameter" (diameter of spread 
of the sonar beam at the surface, a function of beamwidth and cruising 
depth) and the "ping spacing" (horizontal distance between successive 
sound pulses, a function of ping frequency and submarine speed). 
Normally the ping spacing was between 1.3 and 1.5 m, and the surface 
beam diameter over almost the whole track was 3.17 m, implying a very 
narrow beam of less than 3° width. 


Initial processing of the tapes was done at the Arctic Submarine 
Laboratory using a Univac 1108-1110 computer. Corrected depth data 
were merged with positional information to give a final tape in BCD 
format which was forwarded to the AIDJEX Project Office in Seattle for 
analysis. In turn the AIDJEX Office forwarded the tape to the Scott 
Polar Research Institute so that the data could be analysed using the 
same criteria and definitions as those employed in the analysis of data 
from HMS SOVEREIGN (Wadhams 1977a,b; Wadhams and Lowry, 1977). This 
report gives the results of such an analysis. 


2. DATA PROCESSING 


The initial processing at the Arctic Submarine Laboratory 
deleted spurious profile points caused by multiple echoes, fish, air 
bubbles etc., and set these to -10.0 ft. Isolated spurious points 
were then regenerated by linear interpolation. Appendix A describes 
this stage of processing in more detail, including the criteria for 
identifying spurious points. The data file created at San Diego 
consisted of a series of "blocks" each containing about 60 data points. 
The blocks were separated by single lines which usually contained all 
zeroes, but which at periodic intervals contained a position fix 
(i.e. latitude and longitude values). These position fixes thus split 
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FIGURE 1 (above). Route of 
USS GURNARD, 7-10 April 1976. 
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the data file into "intervals" for which the distance travelled is 
known. 


The analysis at Scott Polar Research Institute was carried 
out on the University of Cambridge IBM 370/165 computer using an 
Algol 68C program which is listed and described in Appendix B. The 
program split the input file into "sections". Each section contained 
sufficient intervals to make up 50 km of data, and statistics were 
computed for each of the sections. Figure 2 shows the 27 sections 
involved (with an 18 km end-of-file gap between sections 12 and 13) and 
their precise positions and lengths are given in table I. 


Table I. Positions of the 50 km sections. 


Section Latitude Longitude True length 
kn 
1 Starts 70° 35" 31.2" Dak? 150 Tae 52.67 
Ends 705 49.2 144 13 18.0 
2 Ends Fi Ue ose 144% 13 45.0 50.51 
3 Ends Le om UL. 2 144 14 29.4 DO oe 
uy Ends Fo 20 144 14 58.8 51.01 
5 Ends rips Ry 144 17 52.8 52015 
6 Ends 7 a0 ee U7in 2 144 18 06.0 51.62 
m% Ends ise Las U 144 20 02.4 52.43 
8 Ends of On ey a 144 22 39.6 52.61 
9 Ends a. 0. 144 22 25.8 50.92 
10 Ends 75 08 05.4 144 22 18.0 50.77 
al Ends 75a ora y Py 143 47 29.4 51.46 
1 Ends 74 52 11.4 142" 4321.6 52.358 
13 Starts 74 39 28.8 142 16 50.4 52.40 
Ends 7* #15 43.2 141 20 22.8 
14 Ends 75.92 09.6 140° 25. 20.4 52.04 
15 Ends 73 28 45.0 139 32 24.0 51.63 
16 Ends 73> O4% O4.2 138 43 47.4 52.89 
i 8 Ends (2) 40° 19.2 L358) lS) 7018 51.59 
18 Ends fe e500) 4236 139 49 49.8 52.52 
19 Ends 72 41 28.8 Tat (25) 15956 53.08 
20 Ends 72 42 49.2 142 56 00.6 50.04 
21 Ends (MOG 71 19.2 Ate 6 2615306 50.41 
22 Starts Fae "Fi 22.2 144 48 31.2 52.68 
Ends 72 43 14.4 146 23 51.6 
23 Ends 72 42 48.0 147 59 18.0 52.60 
24 Ends 72 42 05.4 149 32 13.2 51.56 
25 Ends 72 kb (2558 151 02 32.4 50.00 
26 Ends 7 2a e2si Os, 2 152 38 48.6 55-29 
27 Ends 2 (25 eeS2 154 14 32.4 52.98 


The actual depth data from the Arctic Submarine Laboratory came . 
in the form of equally spaced depth points, the spacing being 
unspecified and varying from interval to interval. For every interval 


the Algol program therefore had to calculate an "interpolation length", 
the true spacing between depth points, by dividing the length of the 
interval (calculated assuming a Great Circle track between the position 
fixes for the beginning and end of the interval) by the number of points 
in the interval. The contribution made by each interval to the overall 
statistics for a 50 km section was then always weighted by the inter- 
polation length, so that the resulting statistics are unbiased with 
respect to horizontal length. 


The statistics generated are described and interpreted in the 
remainder of this report, and the full numerical values are given in 
Appendix C. The program produced two types of output file:- 


(a) a line printer listing of the statistics, with details of 
the intervals making up the section (Appendix C); 


(b) a file containing three frequency tables for each section 
(ice draft distribution, level ice and ridges) from which 
histograms were produced using a separate program. These 
appear in the body of the report. 


Figure 3 shows a small part of the depth data plotted out. 
The two most obvious features are a high-frequency noise superimposed 
on the supposedly smooth ice bottom contour, and the occasional 
shallow depth point occurring within the structure of a pressure ridge. 
The noise is undoubtedly a feature of the recording system and, since 
it is random, it does not have a serious effect on probability 
densities of draft, although it will produce an anomalously high number 
of very small "pressure ridges" in the statistics. The shallow depth 
points are probably real features, caused by the very narrow sonar 
beam probing into fissures and crannies within the loose block structure 
of the ridge. The effect of these points may be quite large in 
causing a single pressure ridge to appear as multiple ridges in the 
statistics. 


3. PROBABILITY DENSITY OF ICE DRAFT 
4, beeDelini tion 


The probability density function P(h) of draft h is defined such 
that P(h) dh is the probability that a random point on the ice underside 
has a draft between h and (h + dh). P(h) should really be expressed in 
the form P(h, Xy t) since it is a function of time as well as of position. 
Further, although P can be stochastically defined at a point x, an 
operational definition requires a profile to be taken over a finite length 
scale in order to arrive at an unbiased estimate of P. This length scale 
must be large enough to give a good estimator of P while small enough for 
the distribution not to change significantly within its compass. We have 
chosen 50 km, but some analyses have been done over shorter (17 km) and 
longer (200 km) length scales, where necessary. 


P(h) is related via the mean density of the ice to the thickness 
probability density function g(h) of Thorndike et al (1975). g(h) is 
important as an input parameter to various models of Arctic Ocean ice 
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FIGURE 3. Plot of 800 m of ice depth data from section 1. 


dynamics and thermodynamics. These include the AIDJEX model 

(Coon et: all, 197) , in which an initial ice thickness distribution 
develops by thermodynamic growth and decay and is continuously 
redistributed by pressure ridge building and wind-driven divergence; 
and the viscous-plastic continuum model of Hibler (19774 in which 
a pressure term is parameterised using the mean ice thickness and 
percentage ice cover. The most sensitive part of P(h) is the thin 
ice component, since it has been shown (Badgley, 1966; Maykut, 1976) 
that most of the heat flow from ocean to atmosphere in the Arctic 
occurs through ice of draft less than 1 m; this is also the ice 
component which is most readily available for ridge building. 


3.2. Results 


Figure 4 shows P(h) plotted for all 27 sections of the 50 km 
length scale, using a depth increment of 10 cm; the numerical data 
on which these plots are based are given in Appendix C. The general 
nature of all the plots is similar:- an initial peak, due to thin 
ice in leads and polynyas; a second, broader peak due mainly to 
undeformed first- and multi-year ice; and a tail due to ice in 
ridges and hummocks. There is some variation from section to 
section, especially in the extent of thin ice present. 


To display these variations more clearly P(h) was integrated 
over four depth intervals, which can be loosely defined as "thin 
ice" (0-0.5 m);. "young ice" (0.5 - 2m); "level ice" (2-5 m) and 
"ridged ice" (greater than 5 m). The separation of types is not 
perfect - parts of ridges, for instance, may appear in the "level 
ice" category - but the categories are indicative of changes in the 
nature of the ice cover. The results are given in Table II. The 
intervals were chosen so as to give a direct comparison with the 
data of Wadhams (1977b) from the heavily ridged offshore zone to 
the north of Greenland and Ellesmere Island. Wadhams found that 
the "thin ice peak" in his probability density functions usually 
occurred at less than 0.5 m draft, hence his choice of intervals, 
but the present results (fig. 4) usually show the peak at between 
0.5 and 1 m, presumably because the profiles were done later in the 
winter (April compared to October for SOVEREIGN) so that the ice in 
polynyas is, on average, thicker. Thus we have also added a O-1 m 
category in Table II to include all of the polynya ice. 


The results show a remarkable consistency of ice conditions 
over most of the experimental area. The exceptions are:- 


a) the percentage of thin ice, which varies over a wide range 
(0.4 to 12.3% for the 0-1 m band) and with no apparent 
consistency of trend. The cause is partly statistical - 
thin ice is contained in a limited number of polynyas which 
are distributed non-uniformly along the submarine track - 
and partly real in the sense that thin ice has a transient 
existence and is constantly being destroyed by ridge- 
building so that changes in the wind field during the 3 days 
of the experiment may cause the thin ice to be radically 
redistributed. 
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Table Il. Percentages of ice cover in different ranges of draft 
(a) 50 km sections 

: 0-0.5 0.5-2 2-5 5 0-1 Mean 
Section sh rt m m+ m draft 

1 Pe yak 47 ual 4.5 5.09 

2 0.4 15 60 25 0.8 4.22 

3 O35 18 62 20 vas) Gre: 

4 Fy @ 78 20 0.4 4.14 

D 1.5 fi 71 21 PX: 3.92 

6 O07 abs Wis: 13 Pris! 3.45 

i 0.9 23 63 ES 9.3 3.19 

8 O22 10 78 1 2.0 bey 

9 0.1 Ve ral 1 O25 3.61 

10 a2 11 74 14 one 3.47 

i 2a) 11 1S 14 4.8 3.40 

ye! he 9 We 18 aD 3274 

13 0.3 g 73 L/ 0.8 3.7 

14 0.2 5 78 16 On, Fis je 

15 0.5 9 78 13 0.8 bs | 

16 0.6 9 79 12 1.0 3°38 

17 0.8 8 ae 19 Sal 3.89 

18 1.8 10 66 22 bre) Brey; 

19 3.5 10 2 15 5.2 3.62 

20 PIaW, 8 ye 18 44 3.69 

2] Dial 17 64 16 8.2 by 

22 035 il m5 14 2.6 5.0 

23 Oe M 76 20 Led 414 

24 0) 9) 18 65 16 dr) 3.65 

25 O25 ' 7) 21 Ley, 4.17 

26 0.6 8 63 28 5.0 4.50 

27 0.8 5 68 af ney 4.61 

Mean 0.9 10 Zh 19 3.4 by 
SOVEREIGN 4.8 4 28 63 6.5 (see 
DREADNOUGHT Pdi T 9 55 56) 655) 4 29 

(b) 17 km sections in offshore zone 

Le O52 9 42 48 0.6 5.58 

Le2 SF 14 47 55) 15.3 "AS 

1G 0.0 9 a), 38 Oo 5.18 

aL O25 1s 56 31 0.9 4.51 

PA Pd 0.3 18 7 24 0.6 4,14 

PRS) On 13 66 a 0.8 4.03 

(c) 200 km sections 

A 0.6 IR) 56 29 bry 4.36 
B (Caribou n-s) 0.8 Vik ral Le Det 3.68 
C 0 11 74 14 2.4 3.46 
D 0.6 8 76 16 V0 3.68 
E Lee 9 fe Ly, age) 4.69 
F (Caribou e-w) 1.6 10 72 17 Ae 3.71 
G 0.6 9 68 23 4.9 425 
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(ii) Ridged ice at the southernmost (and, to a lesser extent, the 
westernmost) end of the profile is significantly greater in 
quantity. The percentage of ridged ice is very high in 
section 1] and diminishes to a fairly steady "equilibrium 
value" by section 3. Clearly the first 2 sections represent 
the "offshore province" of Weeks et al (1971), a heavily 
ridged coastal province where the mean onshore tendency of 
ice drift leads to net convergence and ridge building. To 
a lesser extent sections 26-27 mark the outer edges of this 
province further to the west. Figure 1 shows that sections 1 
and 2 occur over the continental slope (the profile itself 
commences at the 100 m depth contour) and that section 27 
ends just as this slope is again approached off Point Barrow. 
The remainder of the sections are samples of what seems a 
very homogeneous ice cover. 


The mean values over all 27 sections are compared in table II 
with the mean results from SOVEREIGN (a 1000 km profile from 81°N 0°W 
to 84°50'N 70°W) and DREADNOUGHT (a 560 km profile from 85° to 90°N 
at 6°E in the ice of the Trans Polar Drift Stream; Williams et al, 
1975). The best agreement is with the DREADNOUGHT data, although the 
GURNARD data show a somewhat lower mean draft which can be ascribed to 
the beamwidth of the DREADNOUGHT's echo sounder. Clearly the ice 
encountered by SOVEREIGN was far more heavily ridged and thicker than 
even the heaviest section of GURNARD. 


The 50 km gauge is short enough to resolve most real variations, 
but to investigate the rapidly changing ice conditions at the beginning 
of the profile (sections 1 and 2) a 17 km gauge was used, i.e. each 
section was split into three. The results, in table II(b), show a 
steady decrease in the percentage of ridged ice as the submarine travels 
north away from the Alaskan coast. Note the isolated value of 13.3% 
for O-l m ice in section 1.2. 


Finally, to obtain very reliable statistics for large tracts of 
the ice cover, the 50 km sections were combined into a 200-km length 
gauge as shown in figure 2, lettered A to G (A is 150 km only). B 
and F are now the appropriate sections for the crossings under Caribou 
camp, and have the length scale recommended in the AIDJEX model and by 
Thorndike et al (1975). The results for B and F agree very closely. 
Again, the 200 km statistics show that the character of the ice is 
essentially constant over most of the track (B to F), with an increase 
in mean draft and percentage of ridged ice at the western end (G) and, 
particularly, the southern end (A). 


3.3 Statistical reliability 


It is clear from table II(b) that there are progressive changes in 
the nature of the ice cover over sections 1 and 2, and that these sections 
(and probably 26 and 27) differ in nature from sections 3-25. The 
question that remains is whether the variations between sections in 3-25 
are statistical artefacts, i.e. due to finite sampling length, or 
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whether they are due to real variations, albeit minor, in the nature 
of the ice cover. Our null hypothesis is that the ice cover over a 
substantial part of the southern Beaufort Sea (the area sampled by 
sections 3-25) is a homogeneous cover in which any given statistical 
parameter tends to the same value everywhere if sampled over a 
sufficient length of track. 


The question of "sufficient length"is a crucial one. For a 
given sampling length, some parameters are estimated more accurately 
than others. For instance, 50 km of track usually contains enough 
ice types to make it a good estimator of mean draft, and most probably 
of percentage ridged ice. ' It may not be a good estimator, however, 
of thin ice percentage since this ice is contained in a small number 
of polynyas which may not happen to fall uniformly within the length 
gauge. A longer gauge may be required for a good estimate of this 
parameter, and also of such parameters as mean polynya spacing or the 
frequency of very deep pressure ridges. 


We can test for the homogeneity of the data using a non-parametric 
run test (Bendat and Piersol, 2o70). The data is divided into n 
sections for which a given statistic S takes values S; (j=l,n). The 
mean value of S is calculated, and each section is classified as (+) or 
(-) according as S; >S or S;< 5S. The number of runs of consecutive 
(+) or (-) classes in the n sections is found and tested for Significance. 
An exceptionally small number of runs signifies a trend or a clustering 
in S; an exceptionally large number signifies a factor tending to cause 
alternation of high and low S values. The 23 50-km sections (3-25) were 
tested in this way, and an additional test was afforded by 30 17-km 
sections corresponding to sections 3-12 (sections 1-12 were analysed at 
17-km gauge for use in table II). The mean values of the statistics 
of table II over 3-25 and the results of the test are as follows 
(* signifies non-significance):- 


Mean and Test - 50 km Test - 17 km 
standard error n=25 n=30 
0-0.5 m percentage cover LO +80 F2 2.5% * 
0.5-2 m i i 10.3 + 0.9 ry ‘a 
2-5 m i " 72 +OiL * * 
over 5m st i" 16.6 + 0.7 * * 
O- m u : Deut 0.7 rs % 
mean draft " " 3.67 + 0.06 * * 


Only one statistic - the 0-0.5 m percentage cover - was rejected by 
this test, at the 25% significance level and for the 50 km sections. 
This implies that the thin ice percentage does not come from a homogeneous 
population, i.e. that there are significant trends or clusterings in this 
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statistic (which can be seen by inspection of table II) suggestive of 
a process acting with a wavelength much greater than 50 km and 
tending to cause large amounts of thin ice over a number of sections 
followed by small amounts. This process must be the wind stress 
field which causes divergence in one zone of the ice cover and 
convergence in another, on a length scale of hundreds of km. 


Otherwise we can accept the hypothesis that the ice draft 
distributions over 1150 km of track in the southern Beaufort Sea 
(3-25) come from a homogeneous ice cover with constant statistical 
properties. The best values to take for the mean draft and percen- 
tages of ice in various depth ranges are given above, together with 
the standard error. It can be seen that this standard error is 
virtually constant for each class (except 0-0.5 m), implying a greater 
fractional error in the estimates of uncommon ice types than in those 
of common types - a result to be expected. 


3.4% Cumulative probability 
The cumulative probability G(h) is defined by 


G(h) = af P(h) dh 


It is used as a major parameter in the AIDJEX model (Coon et al, 1974) 
where it is known as the "ice thickness distribution". We have 
computed G(h) for the two 200 km sections B and F which cross the 
Caribou camp in the south-north and east-west directions. The results 
are shown in figure 5. 


The two distributions differ slightly but are similar in general 
shape. The median depth kM = 0.5) is reached at 3.2 m and the 
graph is plotted as far as G(h) = 0.99, which is reached at 12.2 m. 
G(h) = 0.999 is reached at 16.4 m. 


4, LEVEL ICE 
4.1 Definition 


Williams et al (1975), in their analysis of the DREADNOUGHT data, 
sought a way of determining the preferred thickness or thicknesses of 
undeformed floes, and the percentage of the ice cover occupied by 
floes of this type. By trial and error they decided that the best 
working definition of "level ice" is that the draft point concerned 
should have draft points 4 m to each side of it differing in depth by 
less than 20 cm, i.e. a local gradient of less than 1 in 40 measured 
on an 8 m gauge length across the point. On account of the high- 
frequency structure in the GURNARD profile (fig. 3), we have relaxed 
this definition slightly and we define a level ice point as one whose 
draft differs from a point 10 m away to either side by less than 25 cm, 
i.e. a 1 in 40 gradient in one direction moving away from the point. 


19 


G(h) 


0-5 i 


P(h) 


1.0 


0.5 


20 


Section F (@-w crossing) 


------- Section B (s~n crossing) 


DRAFT m 


FIGURE 5 (above). Cumulative 
probability distribution 
G(h) of ice draft for the 
two 200 km sections which 
straddle the Caribou camp. 


FIGURE 7 (left). Probability 
density function P(h) of 
level ice draft for whole 
submarine track. Bin size 
0.1m. Dotted lines show 
bins between 2 and 4 m draft 
after normalisation (see 
text). 


4.2 Results 


Figure 6 shows probability density functions of level ice draft 
for all 27 sections at 50 km gauge, again using a 10 cm depth increment. 
Any "level ice" occurring beyond about 5 m draft must be on the bottoms 
of ridges and hummocks. The inferences to be drawn from these plats 
concern: - 


(ae) positions of peak or peaks, indicating preferred drafts 
of level ice of varying ages (young, first-year, second- 
year, etc.); 


(ii) mean draft of level ice; 


(iii) percentage of total ice cover occupied by "level ice". 
These results are gathered in table III. 


There are two methods of estimating Cinta One is to plot the 
probability density function of level ice over the whole track length 
(figure 7); the other is to use a ranking technique on the individual 
50 km plots to detect preferred drafts. In table III this has been 
done by listing in order of size every bin with a probability density 
of greater than 0.5 plus the thin ice peak if it exceeds 0.1. 


It is tempting to ascribe the preferred depths of table III to 
ice of varying ages. Clearly the 2.7-2.8 m category is the most 
popular depth (shown also by figure 7) with 2.3-2.5 m and 2.0-2.1 m 
as lesser peaks. Thick level ice of 3.1 depth and above is confined 
to regions in the far west and far south of the profile, i.e. the 
regions of heaviest ridging. However, figure 7, which at first 
sight confirms the hypothesis of a discrete number of preferred drafts, 
in fact required modification. The 10 cm increments from which 
figure 7 is constructed are not all present with equal probability, 
Since the original draft data from San Diego were quoted in integral 
numbers of decifeet (0.03048 m). Thus some of the 10 cm increments 
contain four such decifeet components, while most contain three. In 
figure 7 those increments with four components have been modified (with 
dashed lines) so that their probabilities conform with the increments 
which have only three components. It can be seen that the 2.4-2.5 m 
peak was really an artefact, while the 2.0-2.1 m, 2.7-2.8 m and 
3.1-3.2 m peaks remain as valid peaks of probability. 


Can we say, then, that 2.0-2.1 m, 2.7-2.8 m and 3.1-3.2 m 
represent, say, the drafts reached by first-, second- and multi-year 
undeformed ice? According to present thermodynamic theories of ice 
growth in the Arctic Ocean, we cannot. The results of Maykut and 
Untersteiner (1971) § as modified and tabulated by Thorndike et al 
(1975), show that ice growing from open water at the end of summer 
will reach a thickness of 1.76 m by April 10 of the following year, 
and 2.04, 2.21 and 2.35 m by April 10 of its second, third and fourth 
years of growth. All these values, and the yearly depth increments 
between them, are less than the level ice values that we have found. 
The identification of level ice draft with ice type therefore remains 
an open question. 
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FIGURE 6. Probability 
density function of level 
ice draft for 50-km sections. 
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Table III. Level ice summary 


% Thin Level ice drafts with 
level ice ice peak probability density >0.5 


RA RA A Na AR a SPR TA ERE ne RSE eS 


Section Mean draft 


] 5.81 55.0 0.6 
2 587 43.5 
3 2.96 joe5 10 
4 3.54 DL ee oe 
#5 33kS ipee a paaek 
6 o.7a 58.8 0.8 Deine ct 
é 2.43 59.8 0.8 Salvaeh ace ed asks f 
8 Ok 58.25 1.0 Veeco. 72.0) 255 2.8 
9 2.88 ave) 2.0 Qettesleose es fle 0 os 
10 2ehe 57.0 Deepest. Me. 6° 255° 2,872, 
al 2h? 56.9 029 aepee 20 225 
12 2.94 54.1 DST A250 52), 42 6°25 
15 2,97 45.5 2 HSA OAD IN ae eet 2 
14 eu 58.1 a woe 2. el 2.9 cet 
15 2.96 59.1 ae et 2sO eek 2.0 
16 2.91 59.3 ep ia sk t256 2952.8" 2-0 
L/ 3.06 53.3 ur U 2.7 2.% 
18 5.09 51.0 Uo PsP. Vike! 
mine) 2.00 58.0 0.3 Qe 2. O 20k OS 
20 oi Vi Ieee Cae af ee 
*21 PAS 54.5 0.8 Ps 
22 2.93 HED aera ee) 6D 
23 552 54. Fol 5-2/5 .5 9-5 
24 2.97 ole | 0.8 51 
25 ard 520 a) o MBN) ot wes) 0) 
26 3.61 47.9 10 
27 3.68 AO 4 Sry 4 
Mean ally: DED 
Mean 2.98 56.14+0.7 
(sections 
3-25) 


* Caribou crossings 
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The values of mean level ice draft in table III are not very 
meaningful, because some ice from keel bottoms appears in the 
histograms of figs. 6 and 7 and biases the mean draft upwards. 
However, the figure for percentage of level ice in table III is 
strongly indicative of the real percentage of undeformed ice present, 
although necessarily an overestimate because of keel bottoms. Again 
we see a low figure for sections 1 and 2 with their heavy ridging; a 
fairly steady equilibrium value for sections 3 to 25 (which passes a 
run test for homogeneity) ; and a final decrease of level ice occurrence 
in sections 26 and 27. The thin ice peak, as shown in fig. 7 and 
table III, is concentrated either around 0.2 to 0.4 m or around 
O: 8s toy. lam. Fig. 7 shows 0.8-0.9 m as the most probable category, 
and the distribution of probabilities in the depth range 0-1.5 m in 
fig. 7 is a measure of the relative frequencies of polynyas of 
varying ages. 


5. DISTRIBUTION OF KEEL SPACINGS 


5.1 Independent keels 


The extent of an independent keel is defined using the criterion 
that the troughs on either side of the keel crest (point of maximum 
draft) must descend at least half way towards the local level ice 
surface, in this case defined arbitrarily as a draft of 2.5 m. This 
is analogous to the Rayleigh criterion for resolving spectral lines in 
optics and is identical to that used by Williams et al (1975), Wadhams 
(1976, 1977b) and Weeks et al (1977) for the analysis of submarine and 
aircraft profiles (for airborne laser profiles the identification of 
the "local level ice surface" is much easier). It differs from the 
criterion of Hibler et al (1974), where the troughs must descend a 
fixed distance (61 cm m for surface ridges) from the peak; Hibler (1975) 
has discussed the effect of this difference in definition on the 
resulting distribution. 


5.2 Theory of spacings 


Hibler eral (1972) showed that if ridges occur at random along 
a track the distribution of spacings between ridges is given by 


Dees) ab @ Jphr ay (px) dx (2) 
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where M is the mean number of ridges per unit length of track and 

Lele 3d dx is the probability that a given spacing lies between x and 

(x + dx) in length. Mock et al (1972) tested this relationship for 
surface ridges using aerial photographs, and found good agreement except 
for an excess of ridges at small spacings. On a purely random theory, 
however, we expect a deficit of ridges or keels at small spacings, on 
account of the so-called "ridge shadowing" effect (Wadhams, 1977b). 

This occurs because keels have a finite slope angle so that their crests 
cannot lie closer than a certain minimum distance Xj yjt- Within this 


distance the shallower ridge is not detected and the ridge-picking 

criterion selects only the deeper ridge. Figure 8 illustrates this effect 

for two keels of relief h, h' (h'>h) relative to the level ice bottom, 

each ridge being of triangular cross-section with slope ~. Under these 

circumstances ; 
Sey ade h “cot. (3) 

A theoretical treatment of the modification of (2) by (3) is complex, 

but an approximate solution is given in a paper awaiting publication 

(Lowry and Wadhams, unpubl.). The net effect is that close spacings 

and shallow ridges tend to be lost preferentially from the distributions 

of spacing and draft. 


5.3 Results 


The complete numerical results for keel spacings are given in 
Appendix C. The results are presented for each of the 50 km sections, 
with a spacing increment of 20 m. Distributions are presented for all 
keels deeper than 5 m and for all keels deeper than 9 m. This is because 
Wadhams (1977b) found that a number of deep floe bottoms appeared in the 
draft range 5 - 9 m and that the theoretical keel draft distribution 
function was valid only beyond 9 m. It was felt, therefore, that by 
taking 9 m as a cutoff a more valid distribution of spacings of "real" 
keels could be obtained. The results from Appendix C were added together 
to yield an overall keel spacing distribution for the whole profile; this 
is shown in figure 9. 


Both distributions (> 5 m and >9 m) show general agreement with (2), 
with the expected deficit at small spacings. This deficit shows itself 
only in the spacing range 0 - 40 m, as opposed to the results of Wadhams 
(1977b, fig. 4), where the deficit extends its influence to 120 m. This 
is probably because of the transducer beamwidth in the SOVEREIGN profile, 
which makes ridges seem broader and less steep than they really are 
(Wadhams, 1977c). At large spacings fig. 9 shows a positive deviation 
from (2), which must be due to an additional effect upsetting the purely 
random distribution. I suggest that this effect is simply the presence 
of leads and polynyas, which interpose occasional smooth stretches of ice 
into the otherwise random icefield and thus generate an anomalous number 
of large keel spacings. 


It should be noted that the lines of best fit to the rectilinear 
parts of fig. 9 do not have a gradient of WM For the 5 m cutoff the 
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FIGURE 8. Illustration of the keel shadowing effect for two keels of 
separation x and relief h,h' relative to local level ice draft. 
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FIGURE 9. Distribution of keel spacings over whole submarine track. Bin 
size 20 m. Results are plotted for keels deeper than 5 m and 9 m, and a 
straight line is fitted to the central portion of each curve. 
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gradient is 5.2 ( = 7.3) while for 9m it is 2.9 

( M= irae One expects a gradient of magnitude greater than M 
because keel shadowing implies that the original "population" has 
been reduced before entering the statistics. This is so for 9m 
but not for 5 m, again implying that 9 m is a better cutoff to use 
so as to obtain a population of "pure" keels. 


Appendix C also presents a tabulation of spacings classified 
according to the depth of the deeper keel of the pair which defines 
the spacing, i.e. a tabulation of x against (h' + 2.5) in figure 8. 
This is an attempt to find the best « for use in (3). In Appendix C 
the tabulation is given for all 27 sections and for 5m and 9m 
cutoff. We have taken the results for 9 m cutoff and plotted them 
in bins of 2 m depth increment and 20 m spacing increment. The curves 
(figure 10) show a peak at a spacing which progressively increases with 
depth. In fig. 10 the position of the highest peak in each curve has 
been plotted against the relevant depth (the heavy black dots), and it 
can be seen that the increase with depth is roughly linear. A line of 
best fit has been drawn through these points which, when applied to (3), 


gives a value of 1323" form Of course, these curves show that there 
is no one value of x, otherwise there would be a sharp spacing cutoff 
within which no keel pairs are to be found. Instead, there is a range 


of x, a range which is spread out still further by the fact that the 
keels are not being profiled orthogonally but at various angles of 
encounter. Wadhams (1977c) dealt with this statistical averaging 
problem. Our value of 13.3° is, in any case, an underestimate 
because it refers to the peak of each spacing distribution rather than 
to the spacing at which keel pairs begin to be found. However it is 
indicative of the validity of the "keel shadowing" concept. Again, 
Wadhams (1977c) found by actually measuring the slopes of keels on sonar 
profiles that the slope angle distribution had a peak in the range 
16-20° but, when adjusted to take account of angle of encounter, the 
mean value of & came to 32°. Thus our 13.3° figure does not appear 
unreasonably low. 


6. - DISTRIBUTION OF KEEL DRAFTS 


6.1 Theory 


The theory of keel drafts which has been most extensively tested 
against observation is that of Hibler et al (1972). They used a 
variational calculation which gives the most likely distribution of 
geometrically congruent ridges that will yield a given volume of 
deformed ice. The result is 


P_(h) dh = 2¥h exp (Aho) exp(-\h>) dh (4) 


where P,(h) dh is the probability that the draft lies between 
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SPACING m 


DEPTH m 


FIGURE 10. Spacing distributions for keels of draft greater than 9 n, 
in 2 m draft increments. The spacing corresponding to the maximum of 
each distribution has been plotted as a heavy black dot. This spacing 


increases linearly with keel relief. 
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h and (h + dh); 
h is the mean draft; 


be is a low value cutoff below which keels are not 
included in the statistics; 


») is a parameter which must be derived by iteration from 


exp(-)h?) =h (x1)? erteee hy) (5) 


This has been shown to give a better fit to submarine sonar 


observations than an alternative distribution proposed on empirical 
grounds by Diachok (1975): 


P.(h) dh = 3 exp (-h?/a) dh (6) 


Bae 
with a = 2h/re, 


Again, (4) is modified by the ridge shadowing effect (Lowry and 
Wadhams, unpubl.), but the modification has only a small effect at the 
low-draft end of the distribution where it causes a slight deficit of 
keels; this is a much less drastic modification than that applied to 
the spacing distribution. 


Recently it has been found that surface ridge sails, to which this 
theory was also thought to apply, actually obey a simpler negative 
exponential distribution of form 


P.(h) dh = B exp (-bh) dh (7) 


with B,b as parameters, provided the sails are identified using the 
Rayleigh criterion (Wadhams, 1976; Weeks et al, 1977).  Hibler (1975) 
showed that the same data can be made to fi tes) or (7) depending on 
whether the Hibler (constant trough depth) or Rayleigh ridge-picking 
criterion is used. 


6.2 Results 


Figure 11 shows the distribution of keel drafts for all 27 
sections at 50 km gauge, expressed as keels per 100 km track and using 
a 1 m depth increment. The first bin (2-3 m) actually contains only 
keels from 2.5 m to 3 m in draft, since 2.5 mis the zero datum for the 
Rayleigh criterion. "Keels" of less than 5 m draft can be assumed to 
consist mainly of the bottoms of undulating floes. 


The overall data from all 27 sections were collected and are 


plotted in figure 12 on a semi-log scale for all keels deeper than 
Pyitie The result is an exceptionally good fit to a straight line, not 
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FIGURE 12. Distribution of keel drafts plotted on a semi-log 
scale for data at 1400 km, 400 km and 50 km length gauges. 
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only for the overall data but also for data at 200 km gauge (e.g. 
section F, also shown in fig. 12) and even 50 km gauge (section 1, 
figs?) This shows that the keel draft distribution obeys the 
simpler relationship (7) rather than the Hibler relationship (4). 
This is a most unexpected result, because (7) was hitherto thought 

to be valid only for ridge sails; the SOVEREIGN keel data, analysed 
using the Rayleigh criterion, follow (4) with a high degree of 
exactness as do all other published sonar profiles. Wadhams (1977b) 
suggested that sails may not follow (4) because they contain only a 
small proportion of the mass of a ridge and their shape is determined 
largely by accident. This, cannot explain the present result. We 
must conclude either that keels in the Beaufort Sea have a different 
nature from those in the Eurasian Basin, which is unlikely, or that 
the apparent distribution of independent keel drafts is dependent in 
some way on the type of sensor employed. 


Profiles of ridge sails are always obtained using a laser 
profilometer, which has a pencil beam capable of recording much of the 
fine structure of the sail, including crevices and troughs between the 
blocks (the limitation being the integration time of the laser electronics). 
The sonar employed by GURNARD also had a narrow beamwidth (its precise 
characteristics being classified) and, as shown in fig. 3, it also appears 
capable of recording the fine structure of a keel, probing into clefts . 
and hollows between the submerged blocks. DREADNOUGHT (Williams et al, 
1975) used a sounder with a very wide beam, and SOVEREIGN (Wadhams, 
1977b) had a sounder with a wide beam in the fore-and-aft plane 73 
and a narrow beam in the athwartships plane (5 ). Wide-beam sounders 
smooth out the structure of a keel so that it is always perceived as a 
single wedge (see, for example, the profiles in Wadhams, 1977c), and 
even the application of reconstruction equations (Williams et al, 1975) 
cannot regenerate this fine structure. Now Hibler's theory depends on 
the concept of geometrically congruent ridges, each an entity of .the 
same shape possessing mass and potential energy which depend only on its 
depth. A wide-beam sounder forces keels to approximate to this 
concept by smoothing out any incidental structure that they may possess 
and leaving them as discrete entities. Thus narrow—beam echo sounders 
and laser profilometers produce one type of ice profile with ridge 
height characteristics obeying (7), while wide-beam sounders produce 
another type, obeying (4). A narrow-beam sounder, by splitting many 
ridges into multiple "ridges", sees a greater ridge frequency than a 
wide-beam sounder (e.g. in Wadhams, 1977b, the sail frequency is a 
multiple of the keel frequency for the same ice cover). 


This hypothesis now covers all results except laser profilometer 
data reported by Hibler et al (1974), which still obeyed (4) but which 
were not analysed on the Rayleign criterion. A crucial test of the 
hypothesis would be to smooth the GURNARD profile artificially by 
convolving it with the beam pattern of a wide-beam echo sounder, and 
to observe the effect on the resulting statistics. We hope to report 
on this computer simulation in a later publication. 


It was shown in Wadhams (1977b) that the parameters B,b in (7) 
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can be expressed as simple functions of h and m. 


in the form 


n(h) dh = B exp (-b h) dh 


where n(h) is the number of keels per km track per metre draft 
increment, then 


co 


N= f n(n) an 
hy 
and 
co 
h = Jb n(h). ah 
Oo 
so that 
= -1 
be (h — he) 
and 


to7( 12). 


B = Ab exp(b hy) 


If (7) is rewritten 


(8) 


(9) 


(10) 


(11) 


(12) 


Thus Mand h are the two parameters of the keel draft distribution 
from which the whole shape of the distribution can be deduced using (8) 


sections, using h, = 5m and 9 m so as to be consistent with the 
statistics of Wadhams (1977b). The results are also plotted in 
figure 13. The following conclusions can be drawn:- 


(a) 


(b) 


The major part of the profile (sections 3-25) has a ridging 
distribution which is very homogeneous and which falls within 
narrow limits of variation. These limits are extremely 
narrow for ho = 5 m (5.4-8.2 forNX; 7.2-7.8 m for h) and 
somewhat wider for hy = 9 m, probably because of the smaller 
number of keels involved. Mean and standard deviation for 
these four parameters are given; run tests show that we can 
accept the hypothesis of a homogeneous ice cover with respect 
to ridging intensity. 


Sections 1-2 and 26-27 fall clearly outside the range of 
variation of the other sections (fig. 13), indicating much 
heavier ridging (greater ) in these parts of the track. 

This result agrees with what we have found from the probability 
density functions; at these two extremities of the track there 
are more pressure ridges per unit length, a greater mean keel 
draft, a greater mean ice draft and a greater proportion of 
deformed ice. 


There is no clear positive correlation between and h, 


although the four heavily ridged sections have both a high M 
and a high h. 
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Table IV shows these parameters tabulated for all the 50 km 


Table IV. Pressure ridge frequencies and mean drafts for 50 km sections 


: Draft>5 m Draft >9 m Maximum 
Becbion No. per km Mean draft No. per km Mean draft draft m 
1 14.77 BDz at Let} PTs’ 
2 10.28 8.30 Bh UN 1220 28.83 
3 Sell 7.69 1.81 11.79 21.64 
hs 7.41 7.02 oD) Liev 19.60 
5 7.92 Te 1.90 11.61 22.65 
6 5.68 1d 1320 baad 22.65 
7 5.49 Fa 0.99 12.06 20.54 
8 5.47 7.20 0.91 11.84 20375 
9 Sel7 Lado 1.43 TAS5L 20.09 
10 6.40 (Roel 1.40 11.34 16.86 
1 6.59 Tae 1.24 Ti.25 22009 
12 7.08 7.65 pists ees: 25.05 
(Short Gui7 7.78 1.06 13.10 51512) 
13 6.60 7.09 1.49 11.89 20.97 
14 6.69 7.68 1.58 12.43 24.93 
15 5.46 (poe 0.89 Beg 25..99 
16 6.24 1-59 Lis Tl ei) 18.96 
iy: ia LD 7-79 1.65 Loe 24.84 
18 7.46 7.64 1.64 11.40 18.55 
19 5.50 Do alee LS 26.73 
20 io 4<50 Foe Lines. 20.36 
21 6.98 (ee Leo Ll 3e 19.99 
22 543 ren Le 37 Ll {Soudan be bel Oe 
a) veOd 7.65 L565 12.04 24.48 
24 6.36 7-49 Lowe L250) 22.01 
25 7 woe 7.66 1.78 A erp 22.07 
26 Snn2 8.25 2105 La Je 29.23 
27 S72 8.34 2.60 12.69 29.14 
Mean of main 6.71 7.55 1.39 Jal eo 
group (3-25) +0.20 +0.03 +0.06 +0.07 
Mean of Caribou 
crossings 6.65 7.60 1.42 11.88 
(4-7, 20-23) 
SOVEREIGN 5.68 11.74 3.74 14.19 4541 
DREADNOUGHT 4,20 9.57 2.00 LZ / 30+1 
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(d) For hy = 9 m both the frequency and the mean draft are much lower 
than those found by SOVEREIGN in the very heavily ridged zone off 
north Greenland. Only section 1 exceeds the SOVEREIGN data in 
keel frequency, though not in mean draft. The data compare very 
well with DREADNOUGHT data from the central Eurasian Basin: 
although it may appear from table IV that DREADNOUGHT data have a 
lower M and higher h, the difference can be ascribed to the program 
of Williams et al (1975) which, by applying a harsh version of the 


Rayleigh criterion with sea level as the zero datum, lost many 
shallow keels from the statistics. 


(e) The maximum drafts are surprisingly low. On the whole the deepest 
keel drafts are found in the four anomalous sections (122, 26-27), 
but the deepest draft of all, 31.12 m, in fact occurred in the short 
18 km portion of track that was omitted from the 50 km statistics. 
The general ridging properties of this portion (table IV) are quite 
typical of the sections surrounding it, so that the keel can be seen 
as an isolated event. This is the only keel deeper than 30 m in the 
entire 1400 km of profile, whereas in the SOVEREIGN profile there were 
45 keels deeper than 30 m in 3900 km of track, 39 of them occurring in 
the 1050 km of "offshore zone" north of Greenland (Wadhams, 1977c). 
By a coincidence the deepest keel in the SOVEREIGN profile, 43 m, also 
occurred as an isolated event in an otherwise lightly ridged section 
of ice cover. 


7. LEADS AND POLYNYAS 


The probability density function of ice draft gives the best measure 
of the occurrence and thickness distribution of the thin ice in leads and 
polynyas, and is especially useful for application to heat budget calcula- 
tions. However it is important for a variety of applications in ice 
mechanics, trafficability etc. to know the frequency and width distribution 
of leads encountered by the submarine. Perhaps the most important 
application is to submarine operations themselves - it is desirable to know 
the mean spacing of leads that are large enough to permit a submarine to 
surface. 


A lead was defined as a continuous sequence of depth points in which 
no point exceeds 1 m in draft - thus a polynya broken up by a small floe 
of broken ice counts as two leads. Lead widths were classified in 50 m 
increments and the results calculated for each of the three files making 
up the overall track (File 1 = sections 1-12; File 2 = 13-21; . File 3 = 
22-97). The results are shown in table V. 


On average the aquatic crow has to swim only about 200 m between leads - 
although this figure varies by a factor of nearly 4 between File 2 and 
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Table V (a) Distribution of lead widths : number of leads 
encountered per 100 km of track. 


Lead width m File 1 Filer2 File 3 Overall 
0-50 38 46 
50-100 
100-150 
150-200 
200-250 
250-300 
300-350 
350-400 
400-450 
450-500 
500 
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(b) Mean distance travelled between leads 


Width m ' Distance between encounters 
0 212 m 
50 10.3 km 
100 24.1 km 
150 38.4 km 
200 67.6 km 
500 237 km 
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File 3. However, the lead which he reaches is likely to be very 
narrow, and few leads exceeded 50 m in width. The "exceedence 
table" in V(b) shows that a submarine which requires a 200 m lead 
for a safe surfacing will have to travel 68 km to find one. In 
fact a submarine trying to surface usually investigates every lead 
wider than about 50 m in case the cross-track dimension is 
sufficient to permit surfacing; 50 m leads occurred every 10 km 
in the southern Beaufort Sea during the period of this experiment. 


APPENDIX A. ARCTIC SUBMARINE LABORATORY PROCESSING TECHNIQUE 


The digitized magnetic tapes were unpacked and processed by the 
Arctic Submarine Laboratory on a UNIVAC 1108 - 1110 computer system 
at the Naval Undersea Center. Time marks were removed, and the 
corrected data was then merged with position information and converted 


to BCD format for transmittal to Project AIDJEX. 


After the initial unpacking, spurious profile points resulting from 
multiple echoes, fish, air bubbles, etc., were eliminated by two correction 
programs. Any points that were deleted were set to -10.0. Multiple 
adjacent deleted points were left at this value, but single points were 
filled by linear interpolation. These were flagged by adding 200 to the 


interpolated value. 


The first correction program deleted any points that were greater 
than or equal to 150 feet. At the beginning of a tape, or following 
one or more zeros recorded by the system, the program searched for the 
first valid data point by computing the absolute difference between 
adjacent points (|X, ~ NG agli and requiring that the difference be less 
than 5.0 feet. For example, if the difference between the first two 
points satisfied this criterion, the program assumed that xX) was correct 
and used it as its current X; value. If, however, the difference was 
greater than 5.0, it deleted Xi» made Xy the current X; value, and 
computed a difference (| X, - x, !). It continued in this manner until it 


found a difference that was less than 5.0, at which time it accepted the 
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current X; value as the first valid data point. 


After the program chose a valid XxX. value, it computed the absolute 
difference between adjacent points and compared the difference with a 
threshold number that was set to 20.0 feet. If the Dd, difference 


(|X, - X..,|) was less than the threshold, the program accepted X. 


a+t i+] 
as a valid point and used it as its next Xx, value. However, if the 
9 equal to |x. - Xe ir 
° The program accepted 


difference was greater than 20, it computed a D 


equal to [X, - X 


and, if necessary, a D, is3 
Xiu as a valid point only if both Dy and D. were also greater than the 
threshold value. If Dy was less than 20, Xiu 


was selected as the next Xx; value. If Dy was greater than 20 but 


ee was less, both Xa and X40 were deleted, and X43 was chosen as 


was deleted, and X. 
i+2 


the next X; value. 


Since the system recorded zeros during periods when it didn't 
receive any data, zero words were not used to compute the differences. 
The program required that if Dd, (|X, - Sd igi was greater than the 
threshold, both D, and D, must be greater than 20 or the data were 


2 3 


rejected. In the case that Dd, was greater than 20 and X40 


the program deleted Xia If Dd, and Dy were greater than 20, but 


x. was zero, both X. and X. were deleted. In either case, or if 
i+3 i+] 1+2 


Xi was zero, the program searched for the next valid X; value by using 
the technique described above (third paragraph). 


Was zero, 


In order to maintain continuity between data records, the program 
included the first three data points from the succeeding record with 


the current record when it computed the differences. 


The second program subtracted a surface offset correction from 
each of the profile points. This correction was determined from an 
analysis of both analog and digital data, and is used to shift the mean 
data points in areas of open water to read zero. A few points that were 
already very near zero became negative. The offset was usually one to 


two feet. 


51 


This program also eliminated spurious points that might have 


been missed by the first program by deleting any points during a 


variable time period that exceeded the maximum possible ice depth 


which was determined by examining the analog records. Points 


deleted by this program were not interpolated and appear as -10.0. 


APPENDIX B. SPRI COMPUTER PROGRAM FOR DATA PROCESSING 


The program is written in Algol 68C. It consists of three 


parts:-— 


(a) 


(b) 


(c) 


the main program, which controls the whole computation, 
reading in data, accumulating statistics and outputting 
results; 


a series of statistics procedures for initialising, updating 
and outputting tables, identifying level ice and ridges etc. 
Note that two data-structures are defined, MODE TABLE and 
MODE TWOWAYTABLE, Kach of these structures contains details 
about the classification defining the table, as well as space 
for storing the accumulated frequencies etc.;3 


various miscellaneous "utility" procedures, which may be 
considered as just an extension to the Algol 68C standard 
prelude. 


52 


£(OT*u9%'56-*$ 04d Qez)9271QlR |VSHTEYZ PUY 


oe 


a 
¢s26pyu Qezy £79Aa1 Gey *foud Qey ZIGVL 


Oo SG 


T* 9% $65—-SSebpzyu Qez)971QeyY 954703 F}U} 
Way C9¢v:T} DOT =:S2e6p;4 qezy 34O ves} 


Tt 
4 


I SN1d uo}}xOeSU 


AD YaAS Tf u=epeul uorysezZIIUd 
JSIEP UOFUSZFIUD YON ST yw=FSIP UOFYUSRIUD 19AGY,) FUG 
(2 


$097*56-£12N271 924)F1GC}Y SITE FZEU} 


VW3YN [O9biT] DOT =212AN971 Qe? 40 KeuZ 


v3IY [O09¢%21}3 DOT =:Joud qezy 40 bess 
$=. 4ARUPU LNI 


HOCSIPZSHFERES SPSITEHZ PUSS OH 
6 (U0L}Z99SU *yANOTLIIS ae § FUL TMaus aus {mMoau) ZUpId 


$zno=:}zNopueYS 


(Fuo7r 4Se4% £321 4521) usod yeas 


* SNYL =sP2z4eP 3S NBHL 


(muOSGGe Du =€S2T Je4Us1) . LON 


$(Gp)24P4=s eu} 
2S v4 =:pejyyeys FONE 
ryma 


ussuyymMou 


aw ee 
= 


(WS) aos 


Sos 2PedueZSYP UOYRIES Wor 
H#eeeUudyZoaS Sssad0udereeyH OG 

CwINIde Jd =C7st Jeauty) LIN STIHM 

‘QJQ=ztuoy}z3asu INI 


UPIMIU FSUIYMAUS,ASITISTIVLS XC CI vat PuUyIMau)zuyud 
sfuoy, ySeY *FP7 FSET 6 4=53SIP 18704 WAY 


$ GN3 
[cesty Fury INIYdd $a -— sused uw) FUyUG 
SQODUSDEA(LVELIZISUIVINIEA + COZ([GELIZ}PUIT)NIEA + (LOTsITJIuipNpen=2Guoy 
SEGUUIES(E STZ LE JAUETINIEA + JGOS(LULEG JUITINIEA + (LSSTJausyrwspnyera=27e1 NISBS 
>: GION (6u90744e1 Away AaYH )= Uusod 4°98 DOUd 
$auly avd (S#iT] 
£ GN3 
((9)SO3%(Ze)SOI*R( [RP)SOD 4+ (SE )UISH(TE)ULS )SOdD4e x (1d/I98T) * 09 * 6TO98eT 
‘og t(/yda( pFuoy—cbuoy )=Q Sogisidecze raze S9ST/IAHTFET=ATE Waa 
feewy UL BIUSFSIP 972413 ZFEAUF Sayndwodse *¥# NIOJS 
2 vay (26u07s74e7* FHucys* (4e% TSU )= 99UCFSIP DOU 
HUNK S$GOT= Fury, esuezys}yp uor1;4z93S WAY 


¢¢(szdutRKuns*xuns 
*Ranogesxan0ge #hmoy1a2qexmozyoq INI 
‘Rbdeusexbausy INI CJ 3gN ‘fbe4s ANI [ §J 354 
*RayptxarptAjusxfusAuywexurw (NI ) LONYLS = SIGVLAVMOML JACOW 
*(wns *payyhrem szydu Jy3y ¢s3zdu tsanoqge mo 12q INI 
‘baus Gyan ([] 45yu *Aypejusuyw INT )LONYULS = AVL J3aGOW 


‘jp /T= pPew& uottazyyyud 139A3} 
A(WI)H# *QIOTH FSTP YOPYIFIUD 19h294 TW3Y 
SOO 9G=12A27 wip *OQduf=eFPl4a wip *909G=jsoud wyp INI 
SQUGUQCT= FeJ suaep qoud INI 
Son — FUL UL OM xx ge eee SUM 64 — SOLVS eK KKK HHA ONTULS 


surpzep=suppUeys 
‘(uado 4se OLOO I g=e(TAFUULYD ZYNOpUueZySs *,,yNOZEP,*p~NOZPEP)uado) 
s(uedo 442 O10 | v= (TFUUeYS UppUeZSS,uU}ZeP,ySUPZPEP)UIGO) 
$;nopuez,s=}no SuypuezyS=uy GTS $Fnozepsuyzep BWIA 
NI936 0 
#°**SOILSTAViS XSralve°?# 


53 


Stet st sf tet ot et oF 


ANMPMNORDHDOBNMINORODDN 


¢$aouezsS}P 7e448RU} SN Id @sueyYS}P USO} zZ998 
@(M2)H# £97 PS OUdUSOOOOOT#®WUCZE PP TYEASFZU J= zu} dysejQu}y Way 
¢ (6u07 #2 ey SHu07 YSOYSZOEY ZSOTY )SIUCTS HP= e2ueZsyp 7FeAsezuy W3Y 
$(6u0743e71 pusod yes $6u0y*sze1 TV3RN 381713 
ANYL =2PFe7Z4SRS 
*(6u0y 7SeY ©7071 3501) usod 328 
(FUP PMaUs QUETMISUS,ZFSEZS FE xjj eusod ou 4 
2 2@SNCDSQ S$3}3Z5}32}3S woud PaIpnysxe FeALFIU} SIUS+4+ 44+ 46 2UTMFIU)ZURZId NSHL 
: pezyvexS LON SAIN 
(wpua ye xpf eusod ou » 
w2SNB2PGQ $234S}FHFS WOUJF PSPNTSKS TYEAYSZUY SLUP+t+te+ee* FUPTMIU ZUM NZHL 
atNisjuo 3d =Ce-t 324 rt 31 
s¢(joud wyp=saz}Joudu 
$(epapeesxs joid wWiPw 
6 AOU IKK REKERE  F®UPTMOU)ZUZUG | Joud WIPCaz}Josdu) 
(sy207qQu 42a 7} Joudu s,5 y90RQUSaz} JOUdU wo) xUZUA 
#°°sa7qge} 54235 3}ePdN puUe eazidaud uo SOYRSPZEES JENduos® e# 
£ gO divs Od 
eVV00CN DY =[ Got J9NiA 
se Ij 


(Gv)a4Ppa=rseuzyy $( Fu IMeu)pess 
¢(,,fuepunog yuNnyoseeecee,sauyyMau)zuUzy4d NGHL 
68 we JY =f SOT JeUs yal 
$( Spy )ups=seury {$( eu; TMeuUu)pess 
‘4% SNId S¥207qQU 
+o 
Go 
I4 
((8vo°ce (OQuc—dE dt uudg>4d)) ONNDY =:[ 213s 0udu jay joud | 
foud WiyP=>e7}Loudu) 
‘1 SNId 21} f04du NSHL 
oat-="d GNV 9636-=*d GI 
$((S5 poupspnyea=d INI OG 
Ci.05 
S(euUy{TMaU/)pPePsI OG 
+ OL 
#eyY2D07Q PesttH FATIH 
€(YAZPUPUS WL TIVAYSILNI wt 2@Uy Mou) zUyUA 
£% SNId WAFUPU 
é G=rsyoo1qu *o=r2e7;Joudu ENT $ax2Js04d ENT [JO4d wWIPST J 
#eezVeAsezU1 ssardo0ude *¥# OA 
(4PWPF SPOULZSYP UO}ZSAS KC FPIULESIP uotz. ses 
HO wiNIdae DU =f eitjJe4u}1) LON STIR 


6(G*OFCodaduT—*eu0cGul) =sch sewwns eSpirs =:[TRuewwns eSpyu 
srRuewuns aspyu * (Kuewuns eipr4 NI Cveti 


S(CLOV TIS LE *OTS I )= S4¥d yeyoeds INT (J 


©(QOLEOZ*OT£O07Z£ COGS UFG2EdS Qey) 27QL4NEMOM, SSETEPAPUY 
© INT COZitt*oz:t] DoW =r60eds qey 40 darts 

‘ INI COZ2T} DOT =:679eds gezy 4D Adaus 

© INI COZ:TJ DOT =2628ds Qey 4O xbeuJ 
‘(GO0T*Ozt oz fo7z* COGS OFGIedS Geez) F21Q927MREMOM, SSETEEFIU} 
© INI COZ? t*Ozgzt} DOT =:G9eds gezy 34O — dass 

© INI CLOZ:t}] DOT =:G7eds qey 340 Adaua4s 

‘ INI COZ:1] DOT =:SG2eds Qey 34O xdbeus 

¢go0eds qezy *Gcoeds Qey BJIGVIAVMOML 


AQOanMeNno 


OS o4 4 og ot ost ot 9 
74 od emt ond OF 9 Od ad 


54 


rae ( 
= 


©(6r2eds Qeyz) 21Q27K eEmMOMmYA JZUujyud 
$C 2uyymeus,SQul3W 6 NVHL Y5dd4d30 SS9Gl Yas eu} YMFUS SU) {M]U) ZUyZUG 
e¢¢Ggseds qez) 273qQe3KReEnOMy YUU 
2 (eugwous,SayuLaW SG NVHL 43d3a3q0 S390I14 wf 2U} Mau) pzupud 
$( a> UY (SYURGOP) os 

eSeanqzean K *¢sauzow uy (sHhuzy,seds) sanzea x oo SPUPYMOU)ZUYUA 
S(u-— HIidSG W3N YAd2a0 AB SONIIVdS BOA Ys SOU} WMOU)ZURIA 
© (00Gbt006 ‘<2 uewuwns ashpyu)yKuewuwns jyujzyud 
©€eooce*oostrxuewwuns a6pysu)yKAuewwns Zuysud 
(Awenp*sa2lpysu qez)sz}eZAP 327Q0}% YUuyud 
(e- SSOQIY LNSGNAdAGNIf PU} TM38U)ZUPUd 
¢(Rwanp*zyoanay Qez)STIe2329P 27Qey YUuyud 
$(u—- BDI WIAD Tat eurymauyzuyud 
$(0)=Awwanp INI C J 
®(s3d ye ,sedstjJoud qez)szIeE}FSP 27qQeR% YUrud 
$(e- AIWS0YUd JDIet2uyyMou)zuyud 

€(PAZUUS,=SYEAUSYUL FOUR, 

© we wp e000 Lae eect © uor43e2S) GNNOY 
6,,=2DULYSIP UO}ZVRZIOS, SouUIYMou)ZUIUd 


NSHL 

oO ¢ erueYsIp uoI4ZdeS JI 

$(uorz2asus,NOTLIAS YODA SDILSILVLiSss Puy IM|Us|auUy Mau )zUy ud 
#e8u0r1z99S YwOof $33351}704RS FENAdZYNO® ex¥ 

sedueZS Ip UOIFT9aS SNI1d +S3IP 1EFCF 


Zweez,eAuazur ssado0ud puse *¥# 


Ii 
(c@huewuns vEpr4at sgGy* OQ 06F[ SAPHPr sus { yoSpr4uAhyKuewwns azyepdn 
e(Tkuewmuns sEpri4* CGGbSedS SE sAGpi uu: TyeSpish yKuewwns ayepdn 
H2eds Qey *EKHardedsus{ JoeredsAh *( G29edsusT Jeesedsx) 31QezZfE MOM, azepdn 
Greds qey *[ Gar edsu:f{ JGasedsfh “ff sasedsuz ,jsgeasedsx) ayqeyhRemomy azyepdn 
£(RuUy Guazurfzanay qez4[79PAIVUSTIIFAG71)391Qe}4 FzePAdn 
*(tesaFpiu qezfp saGpirsuu: | jab 
$¢zuUr duvzursfoud qezyff ayyfoudus{, jay} ad)j)o71q9e}4 ezepdn 


OLS ssaBbpyuu *,wo .*}u} daeazyur GNNOY 
.4P w*( JEDIT #22ULZS}EP YEABS|zUY) GNNOY-: 
FU bos 


6 w{YIASTU ‘¢se6pyau @ejyd sad 4S3yp *£3S1Pp YeAUazU sauyyMou)zZUIUd 


©(JSrPK uor~r4azZIsWD Yanay *Szdu uowswazZyYud 7AeNF4 


6 YAN ITU SY FAAT] SL Az,IJOUdU:sT Joyisoud)SzZUIOd YAaAIY AJyzuspy 


(yu GuazyurgS EMU uoiyvazZRyud Yana 7 HPeuH uoyuazyt4Dd 79A24) GNNIDY 


= JJ ypR uorysazyy4ud 19A971 INI 
*(T=2:S}4du uoOrwezyI4d YANG | CGC=S}CduU UuoLywszZyUTD 13N3}1) 
$(2UY AuaezyuySVY_SIP wuoOysszZy sd 19491) GNNOY 
*S3zdu uoyya}IYS 7SAZY LNI 
syeaaqu inI £19A97 INI [179437 WIP. ti 
£(¥ur duazur *GBadsedsu *Gaosedsh * Gaosedsx 
“CoS *[ seGpysuusy, yoOpriskK *( saGpy,uusy, jJe6pysux) sSfFuyoeds dnyas 
#(}uU duayzuri *Gceosedsu ‘*Gasedsh *«Gasedsx 
‘o0s *(C se6pyuusf Jabpy4sNK *( saGpyuusy JaS5pyuxy) sHuyseds dnyos 
*gEa0edsu *Ggesedsu INI 
*gacedsh *Ga0edsx ‘*Gasedsfh «casedsx {NI (SeFpP}4usT j 
*(sa6pi4ustabp ,4sh* aFpyuxt(azyJouduzstjarpsoudysofpyu R}yzuep) 
¢se6pyuu INI Seb Pat eee ed 4x ANI feSpae wyPsTtj 
*Buoz,=sbuo0y %sey £2 eEGTH57Fe1 ¥SL4 


55 


GNA 
(etasezeDd Buyuado,tuayb6sw susado 442 
LIXD (CaePeopaezoxe abpryu wip,ssa)Ssw  2oxs ebpyu wyp 
LIXD (a P2PFFIWXD YAAGY Wi pytsa)yhswu 29x Yansy WIP 
1L1Xa (FSEP VEFOPS ~H=SIUCRSIP YEFOY os 
CDULIMSUS KKK KE KEK KR KR KKK RRR K NNU JO GNant2U} Meu ) Usd ths. 12 
> GO 
#e*eeuUu01zd9aS SSad0ud pusrecry# 
13 
(spesvouby —- uoyzoes yRHUsY O12 Z444+¢¢n0f QUPTMIU )ZUYUd 3ST 

(saBpyyu qezy)e71Qe} YNd}pNO 

$(yaA21 Qez)aFTQey yYNdzyNO 

t(foud qgez)371Qe}4 FNAZNO 

$zyno,zep=7zzyNopueAZs 

#ea7rIS YNdyMo BSyyUMe# 

’ - OS 
( adUueYS}P UO}RTOS 
7OOT#L 33105 26Py4 Qezy 40 b9a4F)) GNNOYA eR Gd ae} 40 be4Js) OC 
sae6py4u qezyz 3O JU OL } YOS 
we Noes, Jo wHOOT wed *ou 0% Sbauf SFP 4 |aSyTeuUISUs# 
¢¢o9e} suap qouds zane} Qezy)sbauJ 27 Q0}% B1H9S 
s(sejf suep qoudtjoud qez)sbaua} 273qQ83 3970985 


eacequruzosors**oHry asode***es NAF @ Wd BIeEzeztse**zoo GOF GNA *°°°*T I -dSVHeeelI-cSVH 


261 
96T 
sé6t 
vet 
£6T 
cot 
T6T 
oé6t 
68st 
eet 
2st 
98T 
Sst 
vet 
est 
2st 
tet 
ost 
62t 
Pa | 
ZLt 
9Lt 
Sét 


56 


$ (GGS2+t4 ) 


Sazyigoud woul 


*c¢ejAsewuwns=wns ‘f 2 JAsewuns=xewh 
> GIDA (Auewwns INI [ J 33y 


$(9x2 72021 wip O109 | 


$¢ 3S 7v4 it pA => 


*( 357V4 Ne 28 ges 


S$} PK u0y 
= QIOA ( 


#49 72 


21 LNI 


‘uaddn 


¢49M07 INI 


(C du+}; JA-([ JA) Sav 
(C du-3JA-(3]f) Sav 3} 


K @dn OL 
6C > Jfuewwns=AKu 


‘cp JAsewwns=-uywhk INI 3594 


’ wos 


NI938 
‘RK INI (J) =Aueuwns a2epdn jO0Nd 


$ GN3Z 
ao 
Ia 
C3} JA=s[ uJ WeAs4 
V2A21 Gdn <u) 
$t SNIld UYU NSHL 
+337 324291 JI 
1 Gdn =>dusy}) 
= 2ybris yonay 
t=<du—} } 
= #421 324391 1908 oa 
R @dn OL }! 4OS 


+464 129491 YO 


GNNOY =:({ se6prysuuyeabpr4ukK sabpruy=zl safpyuujafpiux 


$( 9x2 2e6p34 wip OL109 


$o=:u 
sa AST=PAR *S4du UOPUazIND Yaanaqy=du INI NIOZE 
$SYPAKR UO’uazr4d YSAS} *£S3ydUuU UOtVSzRzZUD 71290397 INI 
aan CJ} a54u *AR INI EC J) =seuUygod yanay A}r.zuepy 9CUd 
©‘ GN3 
ao 
2xueauh=:hA 
© 3 
Id 
If 
$ ANY =se6py4 | AUKTCFR) |]573 
RR=544yU $y =2:26pPr4y | FUCHA) NIHL 
afprisa 4i 
§ BNHYL =sxew geepunof xeweey NOGHL 
FR > yxauARF 47 3973 
I3 
Is 
(RKF=32y | 44>4A) 5573 
13 
RASaY $35 V4 = 226Rips 
[| e6priyux gyn <safpisu) 
‘TT SNild seFpiuu 
geepewsursuod eGbpriseey¥ NAHL 
c/74u > AAR 4% NSAIHL 
ef&pis 41 : 
$ 3S Wd =:xew #eepunojy urweeg NIH 
FR<yxauRR 47 NIHL 
xew 47 


R | O>3™MeURRK) $9OG2S-[ T+! JeXR3 J Oud=HsyxauhA OG 


I O>A") 
$o=:se6piuu 


I —- e43fJ04ud gdh OGL 3} YOS 
£yxauRA ¢OGZ-([ TJer1ppJoswd=sKA Jvay 
$O00000T=:344 WwaynN {$e6pr2yH INI 
$ any. =:xew *© 34967v4 =:26py 4 00S 

NI93€ 


we°esabprisuu pue s6pryuAt eFpyux SOS PUP 


(UOy4azZIUD YH ayRey Aq) saGpyu Zuapuadapuy sayy _Zuapy® ey 


: GIOA (S®6pyau {NI 354 
saSp yuh sa6pjyux INI CC) 35u *973fF04d INI C J)= Se6pyu AJ yquapy DJOUd 


weeesasunpasoud $3}}2S}}ROSRPS OOH 


Uc 


bc ae ae on Et oe Be oe oe oe 


ABANMPNORDADANMNSFNOR MO 


57 


= GIOCA (4% 2ATEVLAVMOML 


) =27Q9e3Remomy ZUTUd 40d 


: GON3 
ea) 


((4 30 PazuByam sydusUGTt*es)O1ze4 GNNDY us} w © T-ZUxAPp+uU pw 


© 4.0} oe SAIP—TUeAIPp+eurwW *,aFUeY uy S3d fo sfez,Ueduad 


wf PUI FMOU) ZU Bud 


© dG {£3(2/30 bests) Sntd Ss OG 
o4U CL TU woss F YOS 


>t - Cr+} sad yeyoads=zu ‘fy syd yeyoads=qtu INI !$)=:8 Way OG 


tT —- $3d 3ye199adS gdn OL ’ w¥O4S 


£(4 30 anoqet,-—anoqge 


*u—-M07EqQ s:aBueu fo yno asnesaq papnirxe seanqyea fo eou 


oft 40 MOT2Q 
sf PULP TMOU)zZUIUd 


2(4 JO P2zuB sem secu GNNOYU fa= Papnyouy szd Jo *ou pazyyfyam “ 
*(+ 3G PatUuUFyom syduty 4O wns )yoyZeUY QNNOY 


Sys ONTLA ueaoy 


uo} 3O S3dus,= pepnysu?: *szd Jo eou 


ox ® 9ULTMOU) ZUIUd 


(AIP 4,30 Sud ISPAYP wf fu 


© yu} wf XOws ,OR 


wf UWS ,wouls pay yJiSSeyYs sue sanyzeEa 


we F9UI>~MaU )ZUIUd 
St-suxAPP+uUlW=xeW INT 


24+ SO AYPH=AIP $4 4O JSussu €y 4O Uyweuyw INT NISS38 
Hee (PITEDS uF 2Q ZOU DALY SAaNYEA bau} SaunsseP)* e¥ 


> GQIOCA (S$34d 1e8;9edSs INI [£ }*2 B7GVL ) 


(L 
S(vGtstuy duszury([y Joh 


u;deap uyw 


$(urs OLG9 | 


=S711}e}EP 34923 yFUIUd DF0Hd 


u=:sSuysedsu 


o<} dl 


oo mgs \V/ ee 


NIS38 


weeuydop uru << yYydap Jo sefFpiu FyuU0E SUapi suode ex. 
#°*(W3) SyUydsp yoaay sYadeaap o} asedsf pue 
‘(w) saBpyu ueem}z29aq sFur1seds 0} voedsx Ss}ase ex 
> GICA (FU d4ezur FYZY 
*sBursedsu INI fey sasedsAh *tadedsx INI [ J] 354 
*u4deap urw Fway feapysN saSpyux iNT [ 3) =sfuyoeds dnyas 3040 


(C2J]Sfe= ww) xew 
6 oo— ueow 
Cayce wf A Oddné,,03 


e®{ ]S*.= wnwyuyw 
a*LvysStu= SAIANYEA * OU 


wf 19MOY*4,2PHuUePY UY SAaNYeA ayy jo 


* GN3Z 


wo ® (C7 IS*CEIS)OHRZC4 GNNOY 


6 FUYY~MOU 
sof PUYYM2U) ZUZUAC 
NI938 


2 QIOA (4eddn *samoy INI *S INI [ J) =Auewuns yuyud JONG 


: ON3 
ao 
Id 


(AR=suywk y uxywASAR) $$ (RR=zsxewh | xeuUuNchR) 
SARK SNtd_ wAs $f SNdd Au NAHE 
uaddn=sf#FK Qny samo .=cAR 41 


- 


$C} JA=AA INI OCG 


omer ke 
nO 


NAHORDAGHNM SHY 


® 
NO 


SOotTCotoCIGantnaaeae 
EF hd cmt td rd Ht ond Ot Od Od wd ed 


oO 


58 


: ON3 


14 
44UB 38M. SNId P274UByem Sydu 
*£ SNId Ss¥du £(z4B,OemMHRAA) SNHId wns 
$2uByemM SNd [CT + (APPS(UYW-AR)) YSILNA J(4 sO 224s) 3S73 
tT SNAId FACQe NZHL 
xeuwcAR 41719 
t SNId *012Q NIAHL 
uyw>AR 4T 
$c} JA=AA INI oa 
F adn OL ? ¥OS 
‘t—-juxAypP+uUyw=xew INT 
*} JO AyP=H=AyP £% JO fJu=Jsu *}% 4O YUYWeUuyW INT 
> pezyuByom szydu=pazyyubrem sydu £3 4H wns=uns W3g 43xy 
ydu=szdu £3 4O SPaoQge=eaoqe *}% 4OQ MOLSQG=MOFOPQ INI ASH NIOZE 
we AQ suoryzngruzpZUOoS BuryzubBsrom *RKR uy sanqyea Fuysn 374Qe2} saz,epdne eg 
GION (448 Lem FAYaH §% ATavi sau *F LENI [ J)= 219%}% 3930 PdN FO0Ud 


wu 


$ GN3 
% 40 AbBeus yvaId $4 30 xbeus 4YVSTD $4 340 beus YVvATD 
$0=:3 430 s¢du=:}3 40 Auwns 
=r} jO xwns=:}3 4O Reanoge=:}4 49 xXanoqe=3} 4O AmoVeq =: AO XMOT2Q 
t‘Rayp=3} 4 Faye Exap=r% AO XAIP 
‘Rfu=2:32 49 Aju ¢xjfu=H:}y 4oO xJu *Auywes2 4O Auyw $xuywe=ty 4O xXUPW NIOAE 
2GIOA (Raypsxayp 
‘Ajusxjushuywexurw INIT ©} BJIGVLAVMOML ASY ) =971Q04hFeEMOM, 2547233 3U DONC 


¢ GN 

4+ 3O tess yvaTdD 

pezyuhrem s4yxdu=:} 40 wns 

=} dO Fa0Qge=3} 34D MO07VEqQ 

ju Surw=3} 4O Uw NIOFE 
3M) =27928}% 2F547F3F1U} DONG 


‘A,P=i+% 30 A 
2 GIGA (Ayptsusuyw INT 


$ GN3Z 
( 0 wf PUL Mau) zUuyud 
ou 4" 
GO  (% *fAye&xu1 (4 3D bAa4s))}ZUUd OG 
+ j3O xJuU OL XY HOS 
(OT *C 83304 AO AKbBosuj)) PUA 
‘(50 *(1T-Ay)*+ 30 Aaip + % 39 Auywyzusd 
(aur ymMou )zupyud- GG 
+ gO F$u or AY wos 
‘(,Ssapyouanbauds RK sanyzyen wf QUYTMU) PUPA 
‘$ dd (% *Cx23(4 AO xdboesus))zU4d OG 
%* 40-xJu OL *Y 30a 
* bau sayoauanbeul Xf PULTMS|U SoU TMaU ) PURI 
oo (¢ ©(T—*%})%% JO XAIP + % 4O xKUyU)ZUId OG 
7 3O xJu OL *} YOS 


Lar ary} SSONTEN XyFAUTMIUS SU YMaU )ZUyUAd 
£( (4 4O s¢duty JO Rurnsyoyze4 GNNOY £48 ue sw a 
*(% 43O $#dusy 4O KWNS)O}XEY GNNOY Sax Ue SU ry 


*} jO Szdus,papnzouy "syd *ou,s aUy Mau 
#4 40 Ranoqety 4O xeaoQge 
*% 40 f¥m0129Q4}% 40 XMOT®FQGSa» 2 YBYY COF Ky 
nt*ubyy OOF xX £MOY CCF KR MOY OCOY X sPAIPNYOX 3, 9UY TMU )ZUZUd NIOZG 


vlt 
e€2t 
cdt 
T2T 
Olt 
69T 
e9T 
291 
99T 
sot 
vot 
e9oT 
cot 
T9T 
o9t 
6ST 
est 
2st 
9ST 
Sst 
7st 
est 


cSt. 


TST 
OST 
6¢T 
Bvt 
2vt 
oVT 
S¢vt 
vel 
Evt 
2vt 
Tvt 
Cvt 
6et 
ect 
Zet 
cet 
SET 
vet 
eet 
cet 
TET 
Cet 
6cet 
Sct 
Zat 
get 
Set 
vet 
ect 
cot 


‘tet 


oct 
6TT 
Btt 
Zut 


59 


je © eon rYsOgOr °° *SHFY 2s0a°***e2 NAF 8 Wd 2e°seector?e?s09 sor GN3° 2 eeeeeTI—-dSWHse8! I—dSVH | 


ete 

ONZ ztz 

(2UIMIU*,,6666u ) yU PUG tt2 

* gO £3132(4% 30 beats) GNNDYH ILNIYd F OAC ote 

3 30 beau} Gaon OL 1 YOS 602 

$( auyymMau) zuyad god 

ssejf suap qoud INIdd 8 202 

$3 40 $3#d¢uU INIYd S&S $+ JO AIP INIUd ¥ 902 

S(L- + JO due} SO ALP + % SO UYW) INTYD % $% 30 UlW LINIdd ¥ NI9O3G S02 
= GQIOA (+ alavi ) =2198% ¥NdyNno D0Yd +02 

€ 02 

¢ OQN3Z zo2d 

Id Toe 

GO 40z2e} x £3304 3D bass) =20 3104 $0,0243) OG 002 

+ §40. duos uaa 66T 

6(% 49 P2xyuByam sydu *« % 30 AYP) Ie} SUEP qoud =1039eJ W334 N3HL S6t 
Oo ¢< + 40 S7zdu JI 26t 

NI93@a 961 


wee,oe} suep qoud, Saw} y Ry:susep AKz~3173QeQOud OF sbau} szPUusAUaTD® e#¥ S6T 
= GIOA (2e$ suep qoud INI *4 37Gvi ) =Sba4s 271Qe% 8309S D0Ud FET 


€6T 

$ GN3Z cét 

Go T6t 

Id O61 

(EL SnAid + 4o Renoge | Axeuchh “21 ytesnid +740 AMOTveR I Ruyw>Afy) 681 
©(4 SNTd # JO xenoge | xxewdxx 34 T SA Id + SO XMOTEQ | Xupyw>oRX) 4S 13 Sst 
I sna + 3O s4du tAA sntid + AG Auns §xx SN Id + JO. xuns zest 

£t SNId CAX*xHI(% BO bass) - 9st 

SL sata CAI SO Abess) $t SATd (xX4}(4 3O xbeus) cet 

4 + (4 SO fares (4 30 Ruyw-AA)) YBILNS = AX vst 

“TL + (% SO “A}P/(4 SO XUYB—-KX*)) YAILING = *% ANI NAHL ‘EST 

Rweew=>AhAR QnNV Ruyw=<hhkK GnNy x¥xMe@W=>xXxX GNY KMUPW=CKXK YI 2et 

$C y JAHAA *¢[( } Jx=xx INI OG tet 

x @dN OL } YOF ost 

‘tT —- + JO Aeu w 4 43D Rapp + KRujyu= Axew 621 

*t - + 30 x*$u we % AO XAIP + KUPWS XXOW INI e2t 

$3 30 Auyw= Aujyw *} 4D XUyW= KuUzPU INI NI938 2Lt 

#°(426uay7 awes pewnsse sfheuse pue x) oy oLt 


© SIOA (2 STAVLAVMOML 3SBU fASx INT (J) =2719083K Cmomy azepdn 50d SLI 


60 


$(jorznopuezs 


$(J=szNROpueys 


a 


-(qse 


$(4o4u2e 0109 5<f uO O>f) 
$uOo SSv -£336 sav =f LNI Oa 
(tiztGeu,py OL t{- AG 6 GdNh wWOXHs3 } wos 


$(e-e=€ 1J6)=6eu 700¢E £$t=:98} INI | $0=:524 INI NI939 


2 INI (f YVHD C J= 2IFA DOUd 


$(usSZst FSsauyymMau)ZUyud s$3yno=szNOpUeYS Synopueys=J BIQIG ) 


>: GIOA (U4 LNI *2S*TS ONTYLS ) = WeABSW 504d 


$(ZS* TS auy mau )zuz,ud $3no=synopueys S$zynopueys=J BJA ) 


=: GIGA (2S*TS ONINLS ) = Bsw 704d 


©(CQb21a<e) ) 
< INE (9° INI )= xewWw DOUd 


: ONS 


(OO00000C IVIO=<e® GNV O=>F) § V=<Q GNV O=>9) =4 WSN NIDFG 


> Way (9°28 Wa3Yy )= 93384 DOUd 


2 (2 


> OO $ W379 3=:€33> CG 101 3 YOGA 
£2? YVHD CST} $S GdN =1 INI ) 
> YVHD C JCS ONTYLS )= DY dO 


ago cG o=:({ff37" OG e Gan 2 D4 F wO4 aa 
e G@gN tf OL 3 YOS ) 
: QIOA (® Wsy (*] 343yu )= uv3aD dO 


se 
~~ 


-( GO v=: !3® OQ #& GdNh OL ’ wos ) 
> GIOA (® Wau € J 435Y¥ Y= YVAD dO 


ago GO o=:C€€*433" OG e&@ San 2 OL §F ywd4 Oa 
e €dh tf OL } was ) 
: GIOA (® .INI [€*] 334 )= YvV3aTID dO 


ee 
~ 


:( Gd =:€ :7" GQ 8% GdA OL $ uosd ) 
=: GEDA (® INI C3 33a ) = YVAWD GAO 


£$(Q4e0=3 8) 
> IW3yH ¢(Q W3uY ee Ws3y 435yu) = 


$(qQ4e=:e) 
z INI (Q INI *@ ANI 334 ) = SNId “dO 
$f= snqd Olud 


‘ ONZ 
qenda 
ago (CIWM=C13xX)=27eNFbS OG aAenbe AJITHM 1OL 3}! YOS 
(AR Q@dN =1)=21eNdbF® FWDDG $x G@dnN =i ANI 
NI92g 
: tong (**x ywvHd [ })= = do 
¢ GN3Z 
> 


$$ GO (€3}]2)P"94 OG 4 Gi } YwOS 
- $$? YWHD CUstj 
NI92ag 

> YVHD C (4 LNI )= 24P4 304d 


6 8 8 eh wt wt et = 1 og 


ANMENORDAROaENMPSNOKRDA 


61 


ee esoHrYug090r®**OHrHN LSOd****ez NO 8 


FNOUZIM (ParJFy¥SANl ZY ys) 


£((1eEAF461Q 00 


so Mat AN $$ (GU LKM—At+e sar 


S(1eA 


Wd cS°8e2°*ct??? 9909 


gor GN5°°*****I] IT—-dSVHe*e4I I-dSVH 


66 

$ GN3S G6 

Go 26 

I4 96 

(LC: Jo)eUuUt4d $$ GnUL =: P2FR4eRS JST S6 

(uo o>) #UP4d NSHL +6 

pexy4eZS 10 GNV ewCe=—}J> AI OG £6 

% oy 1 4wO4d FAs) 

3S Wd =: pezy4uexs 10NG T6é 

(1S7TOA) FUL Gayv= 3 UBVWHD C— ] 06 

#*°*sseo0uez Furpeeay 68 

sjafa}zur ue szpzutiudeex# NIORG Sg 

=: GIOA (1 *18A ENI )= 4#UY4d DCHd 292 

98 

So( SOC ate lele pe 1 Od & @dn ODF 3 tugs G& 
+ -OQLON IG8 CLINI «lL F*9% ENT) = ERs eeo +e 

Es 

£$¢(782)zuU!r dev INIdd ) Zs 

S$ GIGQA, GEYSENI £8 TLN Ie e=! EN see 1s 

oe 

“OG (C3 33)44Ur4d OG F GIN OAL ' yoda )} 62 
: GIGQA (>? HVHD { J) = ENTde come cs 

22 

$t= ILNTod iia G2 

GZ 

$ GN3 $2 

wial: Ez 

S(ecme=s£ IF £ ODS>TRA) Cie 

+ ANTENA ZU Suyf4a),eANSW I QH=o=tA) ib 

SOV )} tdauv =:iiao ek /®)) SPVUNS =" EN TeeoG OL 
(LECEO>IEA) OC T- AG U4 WwONSs ’ YwoOS 69 
VEA— § OOVTEA)=HUEA 1! $2 HVHD CY¥ETJ] NIOJZS eg 
© veo est YEA INI )= FUE dau 908d Zo 

99 

‘ ON S9 

‘ gO (€£33968)2uryud OG F&F GdN OL 2} HOS $9 
S(e- PuNos IF YpP—uoU,* wa)ASw ‘4ouuea eg 

LIXQ .seu 29 

‘(Sau—-=:seu | fauy) | ae) 

: as a) 09 

OtTeoef=3;se}f sse}ulysou=esoau 6S 


62 


G APPENDIX C. FULL NUMERICAL RESULTS 


The numerical results from the foregoing program for 
all 27 sections are lengthy and have therefore been published 
separately as a Supplement to this Report. Readers who have 
not received a copy of this Supplement with their Report may 
obtain one by writing to the authors at Scott Polar Research 


Institute, Cambridge CB2 1ER, Great Britain. 
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